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Abstract 

We present a non-perturbative study of the equation of state in the deconfined phase of 
Yang-Mihs theories in I? = 2 + 1 dimensions. We introduce a holograpliic model, based on 
the improved holographic QCD model, from which we derive a non-trivial relation between 



the order of the deconfinement phase transition and the behavior of the trace of the energy- 
momentum tensor as a function of the temperature T. We compare the theoretical predictions 
of this holographic model with a new set of high-precision numerical results from lattice 
simulations of SU(iV) theories with = 2, 3, 4, 5 and 6 colors. The latter reveal that, 
similarly to the D = 3 + 1 case, the bulk equilibrium thermodynamic quantities (pressure, 
trace of the energy-momentum tensor, energy density and entropy density) exhibit nearly 
perfect proportionality to the number of gluons, and can be successfully compared with the 
holographic predictions in a broad range of temperatures. Finally, we also show that, again 
similarly to the £> = 3 -|- 1 case, the trace of the energy- momentum tensor appears to be 
proportional to in a wide temperature range, starting from approximately 1.2 Tc, where 
Tc denotes the critical deconfinement temperature. 
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2 



1 Introduction and motivation 



The phase diagram of strongly interacting matter has been studied in an extensive experimental 
program since the 1980's. During the last decade, the main heavy ion collision facilities have 
provided convincing evidence for the existence of a new state of matter, which is qualitatively 
different from usual hadronic matter, and appears to behave as a nearly ideal fluid While 
these results confirm the intuitive theoretical expectation that, at high temperatures or densi- 
ties, asymptotic freedom leads to deconfinement, i.e. to the liberation of colored particles from 
hadrons [2], they also reveal that the deconfined plasma is quite different from a gas of nearly 
free quarks and gluons, and should be rather described as a strongly coupled fluid. This makes 
the theoretical description of the equation of state of this system at temperatures close to de- 
confinement particularly challenging: weak-coupling computations [s] have to be pushed to high 
orders, but the convergence of perturbative expansions in thermal gauge theories is generally 
poor, and the evaluation of high-order terms is complicated by the appearence of severe infrared 
(IR) divergences |4|. The latter reveal the mathematically non-trivial structure of perturba- 
tive expansions in finite-temperature QCD (with terms which are non-analytical in a^), and are 
related to the existence of an ultra-soft, chromomagnetic energy scale, which retains an intrin- 
sically non-perturbative nature, and to long-wavelength modes that are strongly coupled at all 
temperatures. 

As a consequence, numerical computations on the lattice are the main tool to derive the pre- 
dictions of QCD at the temperatures probed in experiments, and in the last few years, various 
collaborations have presented results for the QCD equation of state at vanishing chemical poten- 
tial, obtained from simulations including dynamical quarks at or close to the physical point 
At the same time, high-precision lattice results have also been obtained for various equilibrium 
thermodynamic properties in SU(A^) Yang-Mills theories with a large number of colors N, which 
is a particularly interesting limit, for several reasons}^ First of all, the large-A/^ limit at fixed 
't Hooft coupling A = g'^N and fixed number of flavors Nj j9j provides a natural interpreta- 
tion for some non-trivial features of QCD (such as, for instance, the OZI rule [lO]), and leads 
to a topological classification of Feynman diagrams, in which the dominant contributions come 
from planar graphs. This is suggestive of an analogy with similar expansions in closed string 
theory Moreover, for N ^ oo one expects that all correlation functions of gauge-invariant 
operators factorize, and that the functional integral describing a large- gauge theory should be 
dominated by a single "master" gauge field {I2j ; the translational invariance properties of the 



latter are related to the ideas of large- A^ volume independence 13 



As it concerns the phase diagram of QCD-like theories, it is interesting to note that the 



^By contrast, the progress in lattice simulations of the QCD equation of state at finite net baryon density has 
been slower, due to the existence of a severe sign problem [g]. As a consequence, in this case one often obtains 
useful insight from the numerical study of appropriate effective models, see, e.g., ref. [7j. 

^Some of the surprising mathematical properties arising in the large- A'' limit of a generic quantum theory are 
related to the fact that, generally, the large-A limit can be interpreted as a sort of "classical limit". The meaning 
of this statement is made precise in ref. (s], with the definition of an appropriate basis of coherent states and a 
classical Hamiltonian. The construction, however, can be carried out explicitly, and leads to an exact solution, 
only for certain particularly simple models. 
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large- A'^ limit leads to interesting implications at finite density, including, in particular, a possible 
"quarkyonic phase" [14] . Moreover, this limit is also important for applications of the conjec- 
tured correspondence between gauge and string theories [Is] to study the strongly interacting 
plasma [16| , since these computations are done in the infinite- iV limit. 

For these reasons, it is interesting to understand how much the thermal properties of strongly 
interacting gauge theories depend on the number of colors. Recent lattice simulations of large- 



Yang-Mills theories at finite temperature 17,18 have revealed that, in the deconfined phase, the 



bulk thermodynamic observables are essentially independent of N (except for a trivial propor- 
tionality to the number of gluons), and that SU(3) 19 is close to the large- limitj^ This result 



is particularly interesting from a theoretical point of view, since it provides support to analytical 
studies of the QCD plasma based on the approximation of an infinite number of colors, and, 
furthermore, it can shed light onto the effective degrees of freedom relevant for the plasma near 
deconfinement, and/or rule out possible effective models. 

A different perspective on the hot QCD plasma is based on the study of non-Abelian gauge 
theories in a lower-dimensional spacetime. While the case of D = 1 + 1 dimensions is essentially 



trivial 21, in D = 2-1-1 these theories are characterized by rich dynamics 22 , and share 
many qualitative features with their D = 3 + 1 analogues. In particular, at low energies they 
are characterized by a spectrum of color-singlet states with a finite mass-gap and by linear 
confinement, and they exhibit a deconfining transition at a finite temperature Tc. Moreover, 
these theories can also be studied using techniques inspired by the AdS/CFT correspondence: 
comparing the results obtained from first-principle lattice computations with those derived from 
the gauge/string duality can provide a useful test-bed for the application of holographic methods 
to study strongly coupled systems in Z? = 2 -|- 1 dimensions, such as those relevant for condensed 
matter systems at criticality ^23j. Finally, one further motivation to look at the Yang-Mills 
equation of state in 2 -|- 1 dimensions stems from the observation that, in D = 3 -|- 1, the trace 
of the energy-momentum tensor in the deconfined phase appears to be proportional to over a 



broad range of temperatures 18 , 19 , 24 . This behavior seems to be at odds with the expectation 
from perturbative computations, which would rather predict a logarithmic dependence on the 
temperature. In order to understand the physical origin of this characteristic behavior in the 
physical case of D = 3 -|- 1, it is instructive to investigate whether the same phenomenon also 
occurs for a generic SU(A^) Yang-Mills theory in the lower-dimensional setup, given that the 
theories in 3 -|- 1 and in 2 -|- 1 dimensions have both some similarities and some obvious qualitative 
differences. 

For these reasons, in this work we present a systematic study of finite-temperature SU(A'^) 
Yang-Mills theories m D = 2+1 dimensions. Having already discussed the confining phase of these 
theories in a previous work ^25] , in the present article we focus on the equilibrium thermodynamic 
properties in the deconfined phase. In particular, we investigate the strongly coupled regime (close 
to the deconfinement temperature Tc), where physical quantities cannot be reliably caculated via 
weak-coupling expansions, and compare the theoretical results obtained from two different non- 



Similar findings liave been obtained from large-A'^ simulations at zero temperature, see, e.g., ref. 20 and 
references therein. 
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perturbative approaches: holographic computations based on the gauge/string correspondence, 
and numerical simulations in the lattice regularization. 

The structure of this paper is as follows. First, in section[2]we construct a holographic model, 
which is expected to describe the deconfined finite-temperature phase of strongly coupled SU(A^) 
gauge theories in D = 2 + 1. Then, in section |3| we review the basic properties of SU(A^) Yang- 
Mills theories in 2 + 1 spacetime dimensions, and introduce their regularization on a Euclidean 
lattice. Next, in section |4] we present a set of high-precision results from our numerical lattice 
simulations of these theories, for different values of N ranging from 2 to 6. After discussing 
the extrapolation to the thermodynamic and continuum limits, we investigate the dependence of 
the pressure (p), of the trace of the energy- momentum tensor (or interaction measure, denoted 
by A), and of the energy (e) and entropy (s) densities on the number of colors, and compare 
the prediction for A to the holographic model. Section [5] includes a discussion of our findings 
and their implications. The appendix [A] reports the details of the computation of the lattice 
Stefan-Boltzmann limit in D = 2 + 1 and D = 3 + 1 dimensions. 



2 A holographic model 

2.1 Generalities 

The gauge-gravity correspondence ("holography" for short) [15] successfully reproduces most of 
the salient features of large- gauge theories. Although the first examples of holography involved 
supersymmetric and conformal quantum field theories, the correspondence was soon generalized 
to more realistic examples, including theories with linear confinement and no supersymmetry pG] , 
hence in the same class as QCD. However, these models (sometimes referred to as models built in 
a "top-down approach"), that stem from D-brane constructions in type IIA or IIB string theory, 
generally have an infinite number of undesired scalar operators in their spectrum, arising from 
the Kaluza-Klein modes on the internal extra-dimensions of the ten-dimensional parent theory. 

In the meantime, an alternative "bottom- up" holographic approach has been developed [27] : 
it uses minimal ingredients to model the desired features of confining gauge theories on the gravity 
side, in a more direct and "economic" fashion. Our approach here consists of an advanced version 
of the bottom-up construction, that is known as "improved holographic QCD" (IHQCD) |28|29| — 



see ref. |30| for a review, and ref. 31 for similar constructions in the literature. In what follows, 
we first explain and review the setup of IHQCD. 

Holography generally associates the energy dependence of the field theory with a radial di- 
rection r perpendicular to the D dimensions of the Minkowski spacetime in which the gauge 
theory is defined. Therefore the most economic "bottom-up" approach to a D-dimensional QFT 
involves a {D + l)-dimensional gravitational background. In addition, the holographic correspon- 
dence associates a bulk field to each operator that is relevant or marginal in the IR|^ For pure 
Yang-Mills theory, there are two such marginal operators: the energy-momentum tensor T^j^ and 



'^Typically, the holographic description in bottom-up constructions is only reliable in the IR of the field theory, 
because the far UV region of the background suffers from large curvature corrections. However, the reliable region 



turns out to be quite large in llfQCD models 28 . 
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the gluon operator TrF^/A^"^ where A is the dynamically generated energy scale of the theory 
(for D < 4). The bulk fields that are dual to these operators are the metric g^u{r) and the dilaton 
field </>(r). The dependence of these fields on the radial coordinate r corresponds to the renormal- 
ization group scale dependence of the corresponding operators in the field theory. In particular, 
the profile of (j){r) encodes how the (dimensionless) coupling constant gig = (p' j runs with 
the energyj^ On the other hand, in order to make </)(r) run with r in a dynamical gravitational 
setup, one needs to turn on a potential for it. Therefore, the minimal general relativity action of 
the IHQCD is: 

A = M^~^N^ j d^+^x^ [R - ^{d(t)f - V(</.)] + • • • (1) 

(where the ellipsis denotes some boundary counter-terms that should be introduced to render the 
variational problem on geometries with a boundary well-defined; we will not need the explicit 
form of these terms here). The coefficient ^ is an unspecified normalization constanl|^that will not 
play an important role in what follows. In particular, it can be absorbed into by a redefinition. 
Note that the action is proportional to A^^ and to a positive power of Mp, which denotes a 
"reduced" Planck mass thus, in the large- A^ limit of the gauge theory, gravitational interactions 
are suppressed. For convenience, we keep the normalization of the scalar kinetic term unspecified, 
except that we assume > 0. We also assume that the scalar potential has a single AdS minimum, 
that corresponds to the UV limit of the dual field theory: 



V'(0) 

where ^ is the AdS length scale. 



= 0, V{^uv) = ^^^, (2) 

<P=4'vv 



Vacuum solution: The solution to the action in eq. ([T]) that corresponds to the vacuum of the 
dual field theory is of the form: 

ds^ = bl{r) {dr^ + dxl_^ - dt^) , = 0o(r). (3) 

Here the scale factor of the metric 6o(r) is of the AdS form 6o(^) = ^/r only if the dilaton potential 
is constant: V = D[D — More generally, when V is a non-trivial function of 4>, the function 

6o attains the AdS form only in the UV, i.e. for r — )• 0, while it deviates from AdS in the IR 
limit, i.e. for r — )• oo. Its profile can be determined by solving the Einstein's equations, given 
the potential V{(j)). The IR part of the geometry characterizes the confinement properties in the 



vacuum of the dual field theory. In ref. 29 , the various confining asymptotics were classified. In 
particular, we shall be interested in the confining geometries (with gapped and discrete spectrum) 
of the form: 

6o(r) e~^'''^^"+'", a>l for r oo, (4) 



^This is so, because in string theory the field e*^ couples to the operator TrF^/A** ^ in the same way as the 
dimensionless coupling j ts}~^ does. 



This coefficient was set to 4/3 in ref. 28 , as motivated by embedding the theory in non-critical string theory. 
'^The actual Planck mass is the entire expression A/pA^^'''^"^'. 
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where the integration constant A corresponds to the dynamically generated energy scale, and the 
ellipsis denotes sub- leading terms that are typically logarithmic in r [29]. The parameter a is 
directly related to the laige-cj) asymptotics of the dilaton potential in eq. ([T]). Using Einstein's 
equations, it is straightforward to show that eq. Q follows from a potential of the form: 



V((^) ^ const xc/^'^eV 13-1'^ + ... (5) 

where ^ is the normalization factor appearing in front of the kinetic term in eq. ([T]). Therefore, 
a is a parameter of the action, rather than of the particular solution. 

The glueball spectrum of the field theory can be obtained by solving the fluctuation equations 
of the dilaton and the metric, which can generally be expressed as a Schrodinger equation of 
the form —ip"{r) + Vs{r)ip{r) = m?"ijj{r), where ip denotes a generic glueball wave-function that 
corresponds to normalizable fluctuations of the bulk fields]^ 

Here the Schrodinger potential Vs has the following asymptotics: 

Vs(r) oc for r — )• 0; Vs(?") oc r° for r — )• oo. (6) 

Therefore we observe that the Schrodinger potential is bounded both in the ultraviolet (UV) and 
in the IR limits, hence the glueball spectrum is gapped and discrete, if and only if a > 1 p9) . 
We will take a > 1 in the rest of our discussion. In particular, in the limit of large mass, one can 
use the WKB approximation to write down an approximate expression for the glueball spectra: 

ml Cn''-^, for n > 1, (7) 

where C is a constant that depends on the particular model and on the type of glueball. 

Thermodynamics: The temperature is introduced by Wick-rotating the time direction t ^ it 
and compactifying the Euclidean time: r ~ r + 1/T. 

ds^ = bl{r) {dr^ + dxl_^ + dr^) , (j) = 0o(r). (8) 

This geometry corresponds to the thermal ensemble in the confined phase, that we call the thermal 
gas (TG) solution. When evaluated on the TG solution, the action in eq. ([T]) yields the free energy 
of the dual theory. Neglecting the action counter-terms — which are denoted by the ellipsis in 
eq. ([T]) — , the result turns out to be divergent, due to the infinite volume of the asymptotic AdS 
space. Here we shall adopt a particular choice of renormalization, which corresponds to tuning 
the counter-terms, so that the on-shell thermal gas action vanishesj^ 

^The Schrodinger potential is obtained from the background geometry by fluctuating the bulk fields on the 
given background and performing a field redefinition to attain the Schrodinger form. For example, in the case of 
spin-2 glueball s, it is simply determined by the scale factor of the metric as: Vs{r) — {D — 1) log"(6o)/2 + [{D — 



1) log'(6o)/2] 29 . In subsection 6.2 of ref. 29 , it is also shown that the spectrum is bounded from below, as the 



Hamiltonian that appears in the Schrodinger problem can be proven to be positive (semi-)definite. 

^Note that this does not mean that the gas of glueballs in the confined phase has trivial thermodynamics: the 
non-trivial behavior will be encoded in the determinant of bulk fiuctuations around this solution, hence it will be 
suppressed by 1/N^ with respect to the classical saddle-point solution (that, for the TG, is set to zero by our choice 
of the counter-term action). 
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The deconfined phase of the field theory is described by another solution with the same UV 
asymptotics, the black-hole (BH) background: 



ds"^ = b'^{r) [f-^{r)dr^ + dxl_^ + dT^f{r)] , 4> = (t)(r). (9) 

The scale factor 6(r) and the dilaton (j){r) are generally different from their counterparts in the 
thermal gas, eq. ([s]). The blackness function /(r) can be solved in terms of h{r) using Einstein's 
equations as: 

It is a monotonically decreasing function starting as / = 1 at r = and vanishing at r = r/j. 
The latter corresponds to the event horizon of the black-hole, where the 00-component of the BH 
metric in eq. ([o]) vanishes. The temperature of the deconfined gluonic ensemble is given by the 
Hawking temperature of the BH: 

rj.^_fjjA, (11) 

The location of the horizon determines the temperature of the system. The entropy density of 
the ensemble is given by the Bekenstein-Hawking entropy: 

<r,) = ^^-^ = 4vr [Mp6(r,)]^-i , (12) 

where we defined the entropy density dividing the total entropy by the total spatial volume and 



by the square of the number of colors!^ In order to obtain s as a function of T, one has to invert 



the variables using eq. (11). 



The free energy of the deconfined phase is given by the value of the action in eq. ([T]), evaluated 
on the BH background. As we have fixed the counter-terms for the gravity action in eq. ([T]) by 
the requirement that they cancel the on-shell value of the TG action A{TG), the free energy is 
given hy F = A{BH) — A{TG), hence the terms denoted by the ellipsis in eq. ([T]) can be ignored, 
because they cancel in the difference. 

A practical way to determine the thermodynamics of the gluon plasma is as follows. The 



general relativistic system satisfies the first law of thermodynamics (see ref. 32 and references 



therein), therefore one can obtain the free energy density directly by integrating the entropy. 



Here the functions s(r^) and T(r^) are defined in eq. (12) and in eq. (11). The reason for the 
upper bound of integration is as follows. As the location of the horizon r^ tends to infinity, the size 
of the BH becomes smaller and smaller and for r/j — >• oo the BH and TG solutions coincide [32]. 



Therefore, the difference between the actions of the two solutions, hence the free energy should 

^°More precisely, it should be divided by N'^ — 1, but we work in the large- A'^ limit where the difference becomes 
irrelevant. 
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vanish. One can also notice this by comparing the metric functions b{r) and 60 (^) of the BH and 
of the TG solutions in the limit r/j — )• 00, and observe that the difference b{r) — bo{r) vanishes 



exponentially in rh 132 33 , in the entire region < r < r/j. 



Once the free energy is determined as in eq. ^13^, one can change the variable to T instead 



of rh, using eq. (11), and obtain the free energy F{T). Given F{T), all other thermodynamic 



variables follow by standard thermodynamic identities. In particular, one can demonstrate 32 
the existence of a Hawking-Page phase transition at a finite temperature Tc, if the parameter a 
in eq. Q satisfies the condition: a > 1 |^34J. The Hawking-Page transition corresponds to the 
confinement-deconfinement transition in the dual field theory. 

We are particularly interested in the behavior of the interaction measure A as a function of 
T. We define the normalized interaction measure by: 

^(T) = j^=D^ + ^^, (14) 

where we used the standard thermodynamic relation: A = (DF + ST)/Vd-i- The result can 
be immediately obtained, once the black-hole geometry is found, using the formulas given above. 
This is carried out explicitly in the next subsection. 

2.2 Model construction 

In this subsection, we shall discuss a general holographic construction that describes a normalized 
interaction measure A which decays with the temperature as 1/T in D = 2 + 1 dimensions. A 
more ambitious aim would be to construct a holographic model that fits all lattice data presented 
in section [4| This would certainly be possible by engineering the dilaton potential and fixing 
the parameters of the model, which are the parameters in the dilaton potential, the integration 
constants of the equations of motion for 6(r) and 4>{r), and Mp. This aim was successfully 



achieved in the case of -D = 3 + 1 [18,34,|35j. Here, however, we restrict our investigation to the 
general behavior of the interaction measure and leave a more detailed holographic construction 
to future work. 

For this purpose, one can use the following simple, semi-analytic construction. Instead of 
starting from a given dilaton potential and obtaining the metric functions b{r) and /(r) by 
solving the Einstein's equations, one can simply make an Ansatz for the scale factor 6(r), obtain 



/(r) from eq. (10), and derive the thermodynamic functions using the formulas presented in the 
previous subsection. Although this method is not exact but only approximate it turns out to 
provide a good approximation to thermodynamics in a large range of temperatures, as has been 
discussed in ref. . The general conclusion, that we derive in the following, will be independent 
of the details of this approximation. 

At this point, the question is: What Ansatz should one take for the scale factor 6(r)? To 
answer this question, we first recall that 6(r) tends to its TG analogue bQ{r) both in the UV and 



^^In particular, the solution that one finds in this manner does not solve Einstein's equations for a given dilaton 
potential V{(j)), but rather determines the potential. This potential will generally depend on beyond a certain 
value of (j), where the approximation breaks down [361. 
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in the IR hmits, therefore eq. Q leads to: 



6(r) e~('''^)"+-, asr^oo (15) 

£ 

bir) -f - H as r 0, (16) 

r 

where the elhpsis denotes subleading terms. A simple Ansatz that satisfies both limits is 
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b(r) = -e-('^^)", with a > 1. (17) 
r 

As discussed after eq. Q, the parameter a controls the properties of the theory in the IR. In 
particular, a should be larger than 1 for a confining theory. Interestingly, a also controls the 
nature of the confinement-deconfinement phase transition at T^. The transition tends to become 



a continuous one (in particular, a second-order one) as a — )• 1 33,37 . Here we show that the 
parameter a also controls the decay of the interaction measure in the range of temperatures 
between Tc and some intermediate value Tj ^ Tc- The interaction measure is obtained from 
eq. (14) by numerical integration,^^ Fig. [l] shows the interaction measure A (normalized dividing 
by r%2) 

function of the temperature T (in units of the deconfinement temperature Tc), for 
different values of a, in the case of = 2 + 1 spacetime dimensions. We observe that the curve 
becomes steeper with increasing values of a; in section |4j we show a comparison of the curve 
corresponding to a = 3/2 to the numerical results obtained from lattice simulations, see fig. [sj 
One technicality in our present calculation is the value of Tc- In general it is a number of the 
same order as A: 

Tc = coA, (18) 

with Co being some constant depending on the model. In order to determine cq (for a given A), 
we search for a value Tc, such that for rh = rc the free energy difference vanishes: f(rc) = 0. By 
definition, this point corresponds to the phase transition. The transition temperature Tc is then 



obtained by evaluating eq. ( 1 1 ) at r/j = r, 



2.3 General conclusions from holography 

One can understand the dependence of the slope of the reduced interaction measure on the 



parameter a semi-analytically as follows. From eq. (14), one can write: 

rT 
IT, 



^(T) = ^ = D^ + ^ = ^-^ I s{f)df, (19) 



^^Although this is a well-motivated Ansatz, clearly it is not the most general behavior for the scale factor. In 
particular, in the intermediate r region, a less restrictive Ansatz, parametrized by variables in addition to a, would 
allow for a better fit of the lattice data, especially near Tc- However, as emphasized at the beginning of this 
subsection, such a general construction is n ot our primary purpose in this paper. 



^^The upper integration limit in eq. ( 14 1 can be chosen to be a large enough number, so that the change in the 



integral is numerically negligible if the integration limit is pushed to larger values. 
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Trace of the energy-momentum tensor from holography 
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Figure 1: Comparison of the interaction measure A (normalized by and by A^^) that fohows 
from the holographic construction, for different choices of the parameter a. The solid (orange) 
curve is obtained for a = 3/2, the dashed (brown) curve corresponds to a = 2, and finally the 
dotted (maroon) curve is the result for a = 3. 



where, again, we used the standard thermodynamic relation A = {DF + ST) /Vd-i and we 
rewrote the free energy F using the first law as 

F{T) = - I S{f)df (20) 

(where the choice for the lower bound of the integral comes from the requirement that F{Tc) = 
at the transition). Note that, for extremely high temperatures, the expression on the right-hand 
side of eq. (|19|) vanishes, by asymptotic conformality of the theory. Technically, for very large T 



the leading-order term in a 1/T expansion of the integral cancels the first term in eq. (19). The 
contribution that we seek for is the sub- leading term of the integral in 1/T. 

For large enough temperatures T/A ^ 1, one can use the asymptotically conformal result to 
relate T to r^. This follows from the AdS black- hole expression where the blackness function in 
eq. ([9]) becomes: / = 1 — (vh/r)^ . As the background turns into the asymptotically AdS BH for 
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large T, substituting this expression in eq. (11) one finds: 



T 



D 

Airrh'' 



for T > A. 



(21) 



pression from eq. ( 19 ) as: 



Using this expression in eq. (17) and in eq. (12), one immediately obtains an approximate ex- 

(22) 



A(r) = ci 

where we defined the constants 
ci = An I 



-C2 



D 



D 



(23) 



while Co is defined in eq. (18). Evaluating the integral in eq. (22) for large T yields: 



A(r) 



ciC2a 



{D - a) \T, 



l + O 



(24) 



Clearly, this expression is valid only for a ^ D. The case a = D should be treated separately, 
and one finds: 



A(r) w ciC2 Dlog 



-D 



l + O 



-D 



(25) 



The conclusion of this semi-analytic calculation is that, as T increases, the normalized inter- 
action measure falls off as a function of T/Tc, and the precise shape of the fall-off is a power-law 
determined by the parameter a. In particular, this power is not directly related to the spacetime 
dimensionality D, but rather to the nature of the deconfinement transition, which is determined 
by the value of the exponent a in eq. (Q. 

This is in agreement with known results from lattice computations. In particular: 

• In D = 3 -|- 1 dimensions, for pure Yang-Mills theory one expects linear confinement with an 
asymptotically linear glueball spectrum. Then, eq. ([T]) determines a = 2 and from eq. ( 24 ) 
one expects A ~ l/T^ for some range of temperatures above Tc- This has indeed been 



observed in lattice simulations 18,19, 24 



In D = 2 + 1 dimensions, the confinement-deconfinement transition tends towards being 
continuous. In particular, for N = 2 and 3 it is known to be a second-order transition, 
and for = 4 it may be continuous or a weakly first-order one; a general statement that 
one can make is that Yang-Mills theories in 2 + 1 dimensions are more inclined to exhibit 
continuous or weakly first-order transitions than in = 3 + 1 38 -40 . On the other hand. 



as we discussed above, the nature of the transition in holography is determined by the 
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exponent a. For = 2 + 1 one expectf^ it to be smaller than 2. Furthermore, in order to 
have confinement at zero temperature, a should also be larger than 1, see eq. Q. Therefore 
we expect 1 < a < 2 for the Z) = 2 + 1 theories under study. Indeed, as discussed in sec. |4j 
we found that a = 3/2 gives a very good fit to the lattice results. 

One should be cautious with the various approximations that we made in this section. First of 
all, we adopted a semi-analytic approach in determining the holographic background, by setting 
the scale factor of the metric via eq. (17). In the well-studied case of 3 + 1 dimensions, this 

Guided by these results, we made an 



36 



approximation works quite well, as shown in ref. 
educated guess for the scale factor of the metric. Another approximation is to treat the relation 
between T and r/j as in the case of AdS, eq. (21). Strictly speaking, this is only valid for large 
T, where is small enough, so that the asymptotically AdS region sets in. However, numerical 
studies show that this is also a good approximation in a large range of temperatures, ranging 
all the way from the limit of infinite T down to near Tc 35 , 36 . The final approximation is the 



large- limit. Although we believe that the first two approximations can be justified at least in 
a finite range of temperatures, the validity of the latter cannot be assessed by simple arguments, 
but should be checked through case-by-case studies. In the case of 3 + 1 dimensions, the SU(A) 
Yang-Mills lattice data show that, indeed, the equilibrium thermodynamic quantities per gluon 
are essentially independent of A^ flSl, and agree well with holographic calculations [35,41 . We 
hope that this general conclusion also holds for theories in 2 -|- 1 dimensions. This should be 
checked by a thorough study of holographic models in D = 2 -|- 1, that we plan to pursue in the 
future. In the following sections, after introducing the setup to study SU(A) Yang-Mills theories 
in D = 2 + 1 dimensions on the lattice, we compare the lattice results for the equation of state 
with the predictions of the holographic model that we discussed in this section. 



3 Yang-Mills theories in 2 + 1 dimensions 

The continuum formulation of Yang-Mills theories with SU(A) gauge group in D = 2+1 spacetime 
dimensions can be defined via the Euclidean functional integral: 



PAe-^", 5^= /d^x^TVFj^, (26) 



1 

H 

where (which has the dimensions of an energy) is the bare square gauge coupling, F^^ix) 
denotes the non-Abelian field strength tensor, and the functional integration is done over the 
non-Abelian gauge field j4^(x), taking values in the adjoint representation of the algebra of the 
gauge group. Like in D = 3 + 1, SU(A^) Yang-Mills theories are also asymptotically free in 
D = 2 + 1 dimensions; since is dimensionful, perturbative computations for processes at a 
momentum scale k can be organized as series in powers of the ratio g\lk |42j. 

^*In D = 2 + 1 dimensions, one is tempted to set q = 1 for the continuous transitions for N — 2 and 3, by 
the arguments in refs. [33 37 . However, these arguments hold in the large- A'^ limit, and a should receive 



corrections for finite, and small, values of A''. 
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To make the definition in eq. (26) mathematically well-defined at the non-perturbative level, 
one can introduce a gauge-invariant lattice regularization. The common choice is a regularization 
on a cubic, isotropic lattice A of spacing a, which allows one to trade the continuum field A^{x) 
for an at most countable (and finite, if the spacetime is truncated to a finite volume) set of 
matrices U^{x), which are defined on the oriented links joining nearest-neighbor lattice sites. 
The U^{x) matrices represent parallel transporters on the lattice links, and take values in the 
adjoint representation of the gauge group. Their dynamics is defined by: 

3 



-sf 

■' xeAo=l 

where dUa{x) denotes the Haar measure for each Ua{x), and is the Wilson gauge action [43]: 



xeA l<o</3<3 



1- ^Re'IVC/„^(x) 



2N 

with: /3 = ^, (28) 
9oa 



which tends to in the naive (i.e. tree-level) continuum limit a — )• 0, with corrections O(a^), 
and which, being defined in terms of the trace of the plaquette variable: 

Ua^ix) = U^{x)Up{x + aa)Ul{x + a/3)[/J(x), (29) 

is exactly gauge-invariant at all values of the lattice spacing a. 

The expectation value of a physical observable O on the lattice is defined by: 

<0) = ^ / WW'^U^i^) O e-'^ . (30) 

This quantity is a ratio of high-, but finite-dimensional, finite, ordinary group integrals, and 
can be estimated numerically by importance sampling over an ensemble of configurations of link 
matrices. 



Previous lattice computations [38-40 have shown that, similarly to the Z? = 3 -|- 1 case, 
also D = 2 + \ non-Abelian gauge theories exhibit linear confinement at low energy (i.e. the 
interquark potential V{r) grows as ar at large distances r) and a gapped, discrete spectrum of 
color-singlet glueball states. Furthermore, they also undergo a deconfining phase transition at a 
finite temperature Tc, associated with the spontaneous breakdown of a global center symmetry. 

In order to "set the scale" (i.e., to determine the value of the spacing a as a function of /3) 
we used the accurate non-perturbative results reported in ref. [38], which lead to the following 
expression for the temperature (in units of Tc) as a function of /3: 

T /? - 0.22A^2 ^ 0.5 



Nt ■ (0.357A^2 _^ 0.13 - 0.211/Ar2) 



(31) 



The statistical and systematic uncertainties in this scale determination are set by those on the 
value of the critical temperature over the square root of the zero-temperature string tension a in 
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the continuum limit and in the large- limit from ref. 38 . The latter reports Tc/^/o = 0.9026(23) 
in the — )• oo limit of this ratio, with a finite- correction term proportional to A^~^ (and valid 
down to A'^ = 2), whose coefficient is 0.880(43). So our uncertainty on the temperature scale can 
be estimated to be of the order of 1%, and does not have a real impact on our analysis (hence, 
for the sake of clarity, in our plots we do not show the errorbars on the temperature). Another 
potential source of systematic uncertainties is given by the choice of the physical observable to 
set the scale: since all numerical simulations are performed at finite values of the spacing a, 
the determination of the physical scale using different observables can be affected by different 
discretization artifacts; however, the quantitative effect of the induced systematic uncertainty is 
smah, C>(a2). For a comparison with alternative non-perturbative definitions of the scale, see, 
e.g., refs. (39p0p 4]. Finally, note that, in principle, one could also use a perturbative definition of 
the scale; in particular, for the -D = 3 + 1 case high-order perturbative computations in the lattice 
scheme are available in the literature [45j. However, since in the present work we are interested 
in a temperature regime where non-perturbative effects are expected to be non-negligible, our 
determination of the scale is completely non-perturbative. 

In our simulations, we generated the gauge configurations using code implementing a 3 + 1 
combination of local overrelaxation and heat-bath updates on SU(2) subgroups [46] ; for part of 
our simulations, we also used the Chroma suite 1471. In the following, we denote the cardinality 



of our configuration ensembles as n^onf. Having defined the lattice action via eq. (28), the only 
parameters that fix the physical setup of our simulations are the number of colors A^, the Wilson 
gauge action parameter /3, and the sizes of the lattice along the space- like and time-like directions, 
which can be expressed in units of the lattice spacing a as aNg and aNt, respectively. As usual 
in a Euclidean QFT setup (assuming periodic boundary conditions for the bosonic fields), the 
latter quantity is related to the physical temperature via aNt = T^^. 

The parameters of the simulations that we performed for this work are summarized in table [T| 
we chose Ns Nt, which guarantees a good approximation of the thermodynamic limit [48]. In 
fact, finite-volume effects in the T > Tc phase are known to be strongly suppressed, due to the 
screening phenomenon in the deconfined plasma. Note that, as pointed out in ref. [48] (for the 
D = 3 + 1 case), the thermodynamics of a gas of free gluons is sensitive to finite- volume correc- 
tions, which depend on the product of the linear spatial size of the lattice times the temperature. 
In the limit of very high temperatures, such corrections lead to quantifiable corrections to ordi- 
nary thermodynamic relations. However, at the relatively moderate temperatures probed in the 
present lattice simulations, the numerical evidence from all previous studies (both in D = 3 + 1 
and in D = 2 + 1 dimensions) indicates that screening makes finite- volume corrections essentially 
negligible for simulations on lattices with Ng/Nt > 4. In particular, the accurate numerical study 
of SU(3) thermodynamics in = 2 + 1 dimensions presented in ref. [39j provided convincing 
evidence for the strong suppression of finite-volume effects. On the other hand, deviations from 
the thermodynamic limit in the confining phase are exponentially suppressed by the finiteness of 
the mass gap: if niQ denotes the mass of the lightest glueball, then, typically, lattices of linear 
size aNg > 4/mo are such, that systematic effects due to the volume finiteness play an essen- 
tially negligible role in the lattice computation error budget (which is dominated by finite-cutoff 
effects, and by statistical uncertainties due to the finite cardinality of the sampled configuration 
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ensemble). The results at zero temperature are obtained from simulations on cubic lattices of 
volume (aNg)^. In addition to the simulations listed in table [l| for the SU(2) and SU(4) gauge 
groups we also analyzed the configurations corresponding to the Nt = 6 and Nt = 8 ensembles 
(and their respective T = counterparts) taken from ref. i'25|. 

The equation of state of Yang-Mills theories in D = 2 + 1 dimensions can be easily obtained 
from elementary thermodynamic identities. Let Z{T, V) denote the partition function for an 
isotropic system of two-dimensional "volume" V at temperature T; in the thermodynamic limit 
V — )■ oo, the pressure p is related to the free energy F = —TlnZ by pV = —F, and to the trace 
of the energy-momentum tensor A = T'^^ via: 

^ ^Am. ,32) 



r3 dT Vr3 

As discussed above, deviations from the thermodynamic-limit relation pV = —F due to the 
finiteness of the lattice volume can be neglected at the temperatures investigated in this work. 
Finally, the energy and entropy densities (denoted as e and s, respectively) can be obtained from 
e = A + 2p and sT = A + 3p. 

Our determination of the equation of state on the lattice is done according to the "integral 
method" [49| : the trace of the energy- momentum tensor is extracted from differences of {Un)T, 
the expectation value of the average trace of the plaquette at a temperature T: 

A = ^^^{{Uu)t-{Uu)o), (33) 
so that the pressure is obtained by integration over /3: 



3 
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P=^ d^' {{Un)T - {Un)o) , (34) 



'/3o 

starting from a lower integration extremum /3o corresponding to a temperature sufficiently deep in 



the confined phase. We performed the numerical evaluation of the integral in eq. (34) comparing 
the trapezoid rule with the method described by eq. (A. 4) in ref. [50j, which is characterized by 
systematic errors 0{n^'^). Since our scan in /3 values is very fine, the systematic error related to 
the choice of the numerical integration method has a negligible role in the error budget. 



4 Numerical results 

In this section, we present our numerical results for the basic equilibrium thermodynamic proper- 
ties in I? = 2 -|- 1 SU(A^) Yang-Mills theories with = 2, 3, 4, 5 and 6 colors, and compare them 
to the predictions from the holographic model introduced in section [2] By virtue of asymptotic 
freedom, one expects that in the high-temperature limit the thermodynamics of these theories 
reduces to that of a gas of non-interacting gluons, whose equation of state in the continuum reads: 

|3=(^^^-l)f^' (35) 
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N 


X Nt 


n/3 


/3-range 


'^conf at T — 


n^ont at finite T 


2 


48^ 


218 


[7.2,65.0] 


2 X 10^ 


— 




90^ X 6 






— 


5 X 10^ 




64^ 


157 


[10.5,90.0] 


2 X 10^ 


— 




120^ X 8 








5 X 10^ 


3 


48^ X 6 


101 


[16.4, 146.4] 


1 X 10^ 


8 X 10^ 




64^ X 8 


127 


[15.0, 189.6] 


5 X 10^ 


4 X 10^ 


4 


48^ X 6 


261 


[30.0, 246.0] 


2 X 10^ 


1.6 X 10^ 




64^ X 8 


336 


[39.0, 324.9] 


1.5 X 10"^ 


1.2 X 10^ 


5 


48^ X 6 


150 


[43.5,386.5] 


1 X 10^ 


8 X 10^ 




64? X 8 


46 


[60.0,510.0] 


1 X 10^ 


8 X 10^ 


6 


48^ X 6 


129 


[66.0,561.0] 


2.5 X 10^ 


2 X 10^ 



Table 1: Parameters of the new lattice simulations performed for this work: denotes number of 
colors, Nt and Ng are, respectively, the lattice sizes along the time-like and space-like directions 
(in units of the lattice spacing). 72^3 denotes the number of /3-values (i.e. of temperatures) that 
were simulated, in each /3mm < /3 < Pmax interval; the T = and finite-T statistics at each 
/3- value are shown in the last two columns. For > 2, all T = simulations were performed on 
lattices of size {aNg)^. Our analysis also includes part of the data from the simulations reported 



in ref. 25 



where C(3) ~ 1.20205690316 ... is Apery's constant. On the lattice, eq. (35) is affected by cutoff 
corrections: 

|| = (Ar2-l)^i?,(A^,), (36) 
where the correction factor Ri{Nt) can be either estimated numerically or evaluated analytically. 



order by order in an expansion in powers of 51 . For the Wilson action and the integral 



method that we used in this work, the latter computation yields: 

p (N) 1 ^ ^ ^ ^ 227 1 C(7) , 8549 1 C(9) ,n((^/N)^) ^'^7^ 

^'^^'^ = '+ 4iVfC(3)+^iV?C(3)+l2^iVfC(3)+^ ^ ^^^^ 

(see the appendix [A] for details). However, it is important to stress that the cutoff artifacts 
encoded by this correction factor are suppressed at temperatures close to Tc, and hence we do 
not rescale our numerical results by Rj{Nt). 



Since the right-hand side of eq. (35 ) is proportional to the number of gluon degrees of freedom 
(one transverse polarization state for each of the A^^ _ 1 color d.o.f.), in the deconfined phase it 
is natural to normalize the dimensionless ratios p/T^, A/T^, e/T^ and s/T^ obtained in gauge 
theories with a different number of colors, by dividing them by A^^ _ 1, Note, however, that, 
while this is expected to make the results from different groups collapse onto the same curve in 
the limit of very high temperatures, in principle there is no obvious reason to expect the same 
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to be true also at moderate temperatures close to Tc, where the deconfined plasma is far from 
being weakly coupled (and the thermodynamics could perhaps be dominated by different degrees 
of freedom, with unknown scaling properties with N). 

Figure [2] shows our results for the pressure per gluon, in units of T^, as a function oiT/Tc- The 
plot shows the results that we obtained for the various gauge groups, from simulations on lattices 
with Nt = 6 sites in the compactified Euclidean time direction, and for the corresponding space- 
like volumes listed in table [T| As we shall discuss in the following, simulating finite-temperature 
lattices at this value of Nt turns out to be an optimal choice, given that it allows one to reach very 



precise numerical results for the plaquette differences appearing in eq. (34), while keeping the 
systematic effects due to the finiteness of the lattice cutoff under control. The first, striking result 
manifest from fig. [2] is the nearly perfect scaling of p with N'^ — 1: in the deconfined phase, the 
numerical results of the pressure per gluon collapse on the same curve, for all the gauge groups 
that we simulated (denoted by symbols of different colors). This is analogous to what happens 
in 3 + 1 dimensions [17|18] , and can be clearly contrasted to the behavior in the confining phase, 
where, on the contrary, all bulk thermodynamic quantities scale proportionally to 0{N^), i.e. 
are independent of A^, in the large- A limit ^25j. While in principle it is reasonable to expect 
that the equation of state should be qualitatively similar in all of these theories, the remarkable 
quantitative agreement among different gauge groups that our results reveal is completely non- 
trivial. Generally, for SU(A) gauge theories without quark fields, the leading-order finite-A 
corrections with respect to the large- limit are expected to be proportional to A~^, and hence 
could amount to relative deviations of the order of 10% for SU(3), or even larger for SU(2). On the 



contrary, our high-precision lattice results do not reveal any statistically significant evidence 
of such dependence on N. Figure [2] also reveals that the approach to the continuum Stefan- 
Boltzmann limit (which is about 0.1913 . . . , out of the vertical axis range) is relatively slow: at 
temperatures slightly above 7 Tc, the pressure is still approximately 15% off from the Stefan- 
Boltzmann value, indicating that the plasma is still far from an ideal gas of free massless gluons. 



A qualitatively similar feature is also observed for Yang-Mills theories in 3 + 1 dimensions [18 19 
(for which, however, it is important to observe that the energy scale dependence of the physical 
coupling is different). Finally, figure [2] also shows the prediction (solid orange curve) from the 
holographic model discussed in section [2| for a = 3/2. Since the holographic prediction for 



p/[r3(A2 - 1)] is obtained by integration of A/[r'^(A2 - 1)] over T, based on eq. i32n, we fixed 
the integration constant by imposing consistency of the holographic model with the lattice data 
at high temperatures. 
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^^The only difFerences among results corresponding to different gauge groups can be interpreted as statistical 
fluctuations, and/or in terms of the small systematic uncertainty related to the scale setting, as discussed in 
section [3] (for the sake of clarity, the horizontal errorbars associated with the accuracy limits on our temperature 
determination are not displayed in the figures). 

^^Note that the approximations involved in the construction of the simple holographic model considered here 
lead to some deviations from the lattice data for temperatures close to the transition region, and, in particular, to 
an unphysical non- vanishing value for p/[T^{N'^ - 1)] for T T+ . Since the pressure is a continuous function of 
the temperature, this would imply a non-vanishing value for p/[T^{N'^ — 1)] also for T — ^ T~ , in clear contradiction 
with the fact that the number of physical degrees of freedom in the confining phase scales like 0{N^), not like 
O(N^), at large N. 
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Figure 2: The pressure per gluon degree of freedom, in units of T^, as a function of the temperature 
(in units of Tc), for the gauge groups SU(2) (brown), SU(3) (black), SU(4) (green), SU(5) (blue) 
and SU(6) (magenta). The plot shows the results obtained from simulations on lattices with 
Nt = 6. The solid orange curve is the corresponding prediction from the holographic model 
discussed in section [2| for a = 3/2, as obtained by numerical integration of A/[T^{N'^ — 1)] over 
T, with an integration constant fixed by imposing consistency of the holographic model with the 
lattice data at high temperatures. 
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Figure 3: Same as in fig-[2| but for the trace of the energy-momentum tensor per gluon, in units of 
T^. The sohd orange curve is the corresponding prediction from the holographic model discussed 
in section [2| for a = 3/2. Note that, in principle, the holographic model could be refined, to 
match the lattice data also at temperatures lower than 2.5 Tc, through an appropriate choice of 
the dilaton potential. We postpone a more detailed discussion about this issue to future work. 



Figure [3] shows our results for the temperature dependence of the trace of the energy-momentum 
tensor A, normalized in units of and per gluon. The data shown in this plot are the same 
that we used to evaluate the pressure in fig. [2| hence in the deconfined phase they show the 
same, approximately perfect, proportionality to A^^ — 1. For this observable, however, one also 
clearly sees that, in the confining phase (in which the number of physical states is independent of 
— except for the special case N = 2: see ref. |25| for a discussion), the results corresponding to 
different gauge groups do not follow this proportionality law. This effect was not clearly visible 
in fig. [2| because the pressure and the interaction measure are related to each other by eq. ( 32 ) , 
and the signal for p at T < Tc is much smaller than at T > Tc. 

Since p is obtained by numerical integration according to eq. ( 34 ) , whereas A is directly related 
to the plaquette expectation values, see eq. (33), it is most natural to compare the lattice results 
and the predictions of our holographic model for the interaction measure. This is shown by the 
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Figure 4: Same as in fig. [2| but for the energy density per gluon, in units of . 

orange line in fig.[3| which corresponds to the prediction for a = 3/2, as discussed in section[2j As 
one can see, our holographic model accurately captures the non-trivial temperature dependence 
of A for all temperatures T > 2.5 Tc. This good quantitative agreement breaks down in the region 
of temperatures closer to Tc, where the holographic curve falls below the lattice results, and the 
agreement is only qualitative. While in principle the holographic prediction could be adjusted 
to fit the lattice data over an even broader temperature range (by including more terms in the 
dilaton potential), we emphasize that, even for the simple setup discussed here, the model already 
gives a quantitatively correct description for the high-temperature fall-off of A/[T^[N'^ — 1)] and, 
as discussed in section [2j it relates it to the nature of the deconfinement phase transition. 

Our lattice results for the other two equilibrium thermodynamic quantities (the energy and 
the entropy density) are shown in fig. |4] and in fig. [5| respectively. Being linear combinations of 
the pressure and the interaction measure, these quantities obviously exhibit the same, accurate 
proportionality to A^^ — 1 as p and A, and reveal the same type of mismatch between lattice data 
and the holographic prediction for temperatures close to Tc. 

Another interesting problem that we investigated in our simulations is the following: In D = 
3 + 1 dimensions, several authors observed that, in the deconfined phase, the trace of the energy- 
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Entropy density 




Figure 5: Same as in fig. [2j but for tlie entropy density per gluon, in units of T^. 
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momentum tensor appears to be proportional to over a rather broad temperature range 
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19 24|. Since several different interpretations have been proposed for this phenomenon [24], it is 



interesting to investigate whether a similar effect also occurs in D = 2 + 1 dimensions. We address 
this issue in figure [6| by showing our results for the dimensionless ratio A/[T^(A^^ ~ 1)]; plotted 
as a function of Tc/T: if, in the temperature range under consideration, the interaction measure 
is dominated by a contribution proportional to T^, this should result in a linear behavior in the 
plot. This is indeed clearly seen in the figure, hence we confirm that, similarly to the D = 3 + 1 
case, also in D = 2 + 1 dimensions there is a large temperature interval, starting from the value 
where /\/[T^{N'^ - 1)] has its maximum, in which the interaction measure of Yang-Mills theories 
exhibits a quadratic dependence on T. This figure also shows that, at least in the temperature 
range T > 2.5 Tc, the holographic model captures this type of temperature dependence very well. 

The implications of this result are twofold. On the one hand, our finding can be useful to 
shed light on the nature of the phenomenon in D = 3 + 1. In particular, due to the qualitative 
differences of Yang-Mills theories in 2 + 1 versus 3 + 1 dimensions, our result might help to 
rule out some mechanisms that have been proposed to explain the phenomenon in D = 3 + 1, 
if they are expected to be at work also in the lower-dimensional case. On the other hand, 
our holographic model leads quite naturally to a power- law decay of A/[T^{N'^ - 1)] with the 
temperature, and, even more interestingly, it suggests a connection between the order of the 
deconfining phase transition, and the exponent of such power-law decay. As discussed above, the 
fact that, in general, the deconfinement transition tends towards being more discontinuous when 
the spacetime dimensionality increases from 2 + 1 to 3 + 1, can thus be directly related to the 
change from a l/T to a l/T^ fall-off for A/[r^(iV2 - 1)]. 

Finally, we conclude this section with a discussion of the finite-cutoff effects affecting our 
lattice results. As we mentioned, most of the results presented in this paper are based on finite- 
temperature simulations using lattices with Nt = 6 sites in the compactified Euclidean time 
direction. One may wonder, whether the corresponding results are close enough to the continuum 
limit or not. According to our analytical expansion of the Ri{Nt) factor in eq. (37), it turns out 
that, for this value of A''^, the lattice Stefan-Boltzmann limit (evaluated with the integral method 
and the Wilson action on an isotropic cubic lattice) differs from the value in the continuum 
by approximately 5%: an effect much larger than the statistical uncertainties and the other 
systematic errors affecting our data. Thus, in principle one may be tempted to rescale all our 
lattice results by dividing by Ri{Nt). Since Rj does not depend on A^, this would not change the 
fact that the thermodynamic quantities are nearly perfectly proportional to the number of gluons, 
but would lead to slightly different (smaller) numerical values for p, A, e and s. However, in the 
temperature region investigated in this study, this naive rescaling of the results would not be 
correct: the physical reason is that the distortion of the Stefan-Boltzmann limit encoded by Rj is 
due to modes near the lattice cutoff, and those are not relevant for the physics at temperatures of 
the order of Tc. For this reason, we chose not to rescale our numerical results by Rj, but rather to 
repeat our simulations (except for the computationally most demanding gauge group SU(6)) at 
the same temperatures and at the same space-like volumes, on finer lattices, with Nt = 8. Given 
that the leading discretization effects of the Wilson action are 0{a'^), this corresponds to reducing 
the lattice artifacts by approximately a factor 2. From eq. (34) (in which the plaquette mean 
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Figure 6: Similarly to what happens in D = 3 + 1 dimensions flS], in the deconfined phase 
there exists a temperature regime, in which the trace of the energy-momentum tensor A appears 
to be proportional to T^. This is exhibited very clearly by the linear behavior of the data 
displayed in this plot, showing the dimensionless ratio A/[{N'^ — l)r^], as a function of Tc/T, for 
temperatures (approximately) starting from 1.1 Tc (near the maximum in figure [3|. The figure 
shows the results of our simulations on lattices with Nt = 6, with the same color code as in fig. [2j 
and the corresponding prediction from the holographic model (solid orange line). 
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values are always 0(1), for any Nt) it is also easy to see that, when Nt is increased, the difference 
appearing in the integrand on the right-hand side is affected by a fast decay of the signal-to-noise 
ratio. As a consequence, it would become extremely difficult to get sufficiently precise results 
from much finer lattices. Fortunately, however, the discrepancy between our Nt = 6 and Nt = 8 
results for the equilibrium thermodynamic quantities considered in this work appears to be very 
small, as figure [7] shows. This holds in the whole range of temperatures that we studied, and for 
all values of N from 2 to 5. Since our Nt = 6 and Nt = 8 data sets yield compatible results, we 
did not attempt a continuum extrapolation of the thermodynamic quantities, and we can safely 
state that, to the level of precision we reached, our results from the A^^ = 6 lattices are already 
compatible with the continuum limit. 



5 Conclusions 

In this work, we presented a non-perturbative study of the equilibrium thermodynamic properties 
in the deconfined phase of (non-supersymmetric) SU(A^) Yang- Mills theories in 2 + 1 dimensions, 
using holographic computations and numerical simulations based on the lattice regularization. 
This allowed us to combine the advantages of both tools: the former enables one to gain analytical 
insight on the dynamical properties of these strongly coupled systems, while the latter (once the 
thermodynamic and continuum limits are taken) provides numerical results obtained from an ab 
initio approach, directly based on the microscopic definition of the theories for any number of 
colors, without any assumption or uncontrolled systematic uncertainty. 



First, we introduced a holographic bottom-up model, inspired by the IHQCD model 28 



30 , which describes the non-trivial dynamics of these strongly interacting non-Abelian gauge 
theories in the large- A^ limit. This model reveals a non-trivial relationship between the order of 
the deconfinement phase transition, and the dependence of the trace of the energy-momentum 
tensor A on the temperature. In particular, for non-Abelian gauge theories in 2 + 1 dimensions 
(which, typically, are characterized by a tendency towards a second-order or a weaker first-order 
transition than in 3 + 1 dimensions), at temperatures of the order of Tc the model favors a behavior 
approximately compatible with a 1/T decay for the dimensionless ratio A/T'^. 

Then, we defined the non-perturbative regularization of SU(A^) Yang-Mills theories on a 
(2 + l)-dimensional Euclidean lattice, and performed a set of high-precision numerical simulations 
to study their equation of state at T > Tc. We compared the results obtained for different 
numbers of colors, up to A^ = 6, and found that the trace of the energy-momentum tensor and 
the related bulk thermodynamic quantities per gluon are independent of A^, reflecting a strikingly 
accurate scaling of the equation of state in the deconfined phase, over the whole temperature 
range that we probed (up to about 7.5 Tc). This holds for all the gauge groups that we studied, 
including SU(2), and — at least for these equilibrium thermodynamic observables — supports the 
potential quantitative relevance of analytical computations relying on the large- A^ limit (including, 
in particular, those based on the gauge/gravity correspondence). We also found that, in all 
the theories that we simulated, A exhibits a clear, characteristic quadratic dependence on the 
temperature. Both these findings are analogous to those which have been obtained for non- 
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SU(2) 



SU(3) 




Figure 7: Cutoff dependence of the trace of the energy- momentum tensor in units of T^: the 
plots show the results obtained at temperatures ranging from approximately 0.7 Tc to 7.5 Tc, 
from simulations on lattices with Nf = 6 (circles) and 8 (triangles), for four of the gauge groups 
considered in this work: SU(2) (top left panel), SU(3) (top right panel), SU(4) (bottom left panel) 
and SU(5) (bottom right panel). 
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Abelian gauge theories in 3 + 1 dimensions [18, 19 ,24 . Finally, we compared the lattice results 



with the prediction of our holographic model, finding good quantitative agreement, at least for 
temperatures not too close to Tc. 

In the future, we plan to extend the present study, by investigating the holographic model in 
more detail, and by looking at different observables, which could be compared with the results of 
lattice simulations. 
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A Lattice cutoff corrections to the Stefan-Boltzmann limit in 
2 + 1 and 3+1 dimensions 

In this appendix, following the calculation in ref. [Sl], we derive the correction to the Stefan- 
Boltzmann limit due to cutoff effects on the lattice, for the Wilson discretization of SU(A^) 
Yang-Mills theory in D = d+1 dimensions, for d = 2 and 3. Our goal is to evaluate the first few 
terms of the correction to the Stefan-Boltzmann limit, in an expansion in powers of , where 
Nt denotes the number of lattice points in the Euclidean time direction. We take Ng, the number 
of lattice sites along the space- like directions, to be infinite (corresponding to the thermodynamic 
limit). 

Throughout this appendix, we work in lattice units, i.e., we set the lattice spacing a to unity, 
and denote the pressure as p, the spatial volume as V , and the temperature as T. Moreover, in 
the following, we use the = notation to mean equality of two quantities, up to terms which are 
negligible to the order of precision of our computation. 

Notation used throughout this calculation includes: 



Ntsin{pi/2), 
so that: 

Q ip dq y^ 3u^ hy^ 

~ 9Nf ~ lOiVf ' ^ 27iVf ' ay~~iVf 4iV^ ~ 87V|' 
Elementary identities used in this calculation include: 



A sin^(pi/2), X = 2arsinh(Li;), yi 



y = Ntuj, t = 2y, g 



y ^ W 



m? 20Nf 567Vf ' 
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and arsinh(x) = 
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as well as: 
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Furthermore, we also use the following finite-sum formula [52]: 
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Finally, in the following we use Z^^^^ to denote the partition function for one bosonic, massless 
degree of freedom on the lattice. 

We concentrate on the integral method for the lattice determination of the pressure j49j, in 
which p{T), the pressure at a given temperature T, is defined with respect to its value at T = 0. 
The thermodynamic definition of the pressure reads: 
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and, for an isotropic system, in the thermodynamic limit it reduces to: 



T 

P=ylnZ. 



A.l D = 3 + l 

The pressure can be written as: 
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(A.l) 



Thus, taking into account that in the integral method the pressure is defined w.r.t. to its value 
at T = (obtained as the Nt ^ oo limit): 
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Changing variables to yi = Nt sin(pj/2), and expanding pi as: 



Pi = 2 arcsin 



2^ + 



+ 



Nt 3 \NtJ 20 \NtJ 56 \Nt 



+ 



so that: 
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and therefore (by rotational symmetry): 
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Thus the dimensionless ratio p/T^ can be written as: 
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The last expression can be readily evaluated, using the formulas listed above and integration by 
parts. In particular, the first two terms are: 



2C(4) + ^C(6) + ^C(8) + ^C(IO) 



8Nf 



and: 
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Next, note that, introducing the following polar parametrization for the yi coordinates: 

yi = y sin 9 cos (p 
2/2 = 2/ sin 6 sin 
2/3 = 2/ cos 6 



and denoting c = cos 9, one gets: 
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The remaining terms evaluate to: 
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Plugging these results into eq. (A. 3), one eventually ends up with: 
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eq. (A.4) can be rewritten as: 
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which reproduces the expression given in ref. [51] , and extends it to the next order in powers of 
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A.2 D = 2 + l 



An analogous calculation can be done for a system in Z) = 2 + 1 dimensions, for which eq. (A.l) 
gets replaced by: 
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while eq. (A.2) gets replaced by: 
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Accordingly, the Stefan-Boltzmann limit for the dimensionless ratio p/T^, as evaluated on a finite 
lattice using the integral method, reads: 
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Similarly to the D = 3 + 1 case, each term appearing in the last expression can be evaluated 
separately. In particular, one easily finds that: 
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Next, introducing the following polar parametrization for the coordinates: 
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it is easy to prove that: 
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In contrast to the D = 3 + 1 case, the latter expression cannot be rewritten as a simple power 
series in (vr/A^f)^ with rational coefficients, because, for odd values of x, the ir^/(^{x) ratio is not 
just an integer (or a rational) number. However, note that, defining: 
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oo 



k=l ^ ' 



it is possible to write: 
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